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.^ , Abstract 

^^ I A well known result of P. Flandrin states that a Gaussian uniquely maximizes the 

integral of the Wigner distribution over every centered disc in the phase plane. While 
there is no difficulty in generalizing this result to higher-dimensional poly-discs, the 
generalization to balls is less obvious. In this note we provide such a generalization. 



2 ' 1 Introduction 
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The Wigner quasi-probability distribution was introduced by Wigner [16] in 1932 in order 

K^ ' to study quantum corrections to classical statistical mechanics. Nowadays it lies at the 

^ir ■ core of the phase-space formulation of quantum mechanics (Weyl correspondence) , and has 

t^ ■ a variety of applications in statistical mechanics, quantum optics, and signal analysis, to 

name a few. In this note we consider the localization problem of the n-particle Wigner 

distribution in the 2n-dimensional phase space. We state our results precisely in Theorem 

^^ [ 1 below. 

Equip the classical phase space M?^ with coordinates (x, y) with x,y £ M". The Wigner 
quasi-probability distribution on M^", associated with a wave function t/j G L^(]R"') and its 



complex conjugate ip* , is defined by 

W^x, y) = (2^)-" I ^{x + T/2)i)*{x - T/2)e~'^-y dr (1.1) 

The function W^ possesses many of the properties of a phase space probability distri- 
bution (see e.g., [1]); in particular, it is real. However, W^ is not a genuine probability 
distribution as it can assume negative values. 

The localization problem, i.e., estimating the integral of the Wigner distribution over 
a subregion of the phase space, and the closely related problem of the optimal simultane- 
ous concentration of "0 and its Fourier transform ^0, have received much attention in the 
literature both in quantum mechanics, mathematical time- frequency analysis, and signal 
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processing (see e.g. [H [H O [H O [H [H [TOl [121 El [H], and the references within). Bounds 
on the L^ norms were found in [?]• More precisely, the problem of interest for us is: 



The Wigner Distribution Localization Problem: given a measurable set D C 
find the best possible bounds to the localization function 



t>2n 



£{D):=sup W^dxdy, (1.2) 

^ J D 

where the supremum is taken over all the functions ip G L^(]R"') with \\ip\\2 = 1- 

The quantity £{D) is invariant under translations in the phase space, and under the 
action of the group of linear symplectic transformations (see e.g. J15j). There is no upper 
bound on £{D); it can be infinite. Indeed, there is a ■0 £ L^(M) such that J \W^\dxdy = oo 
[H sect. 4.6]. An example is il^{x) = 1 if — ^ < x < ^ and V'(a;) = otherwise. On the other 
hand, the L^ norm of W^ is bounded [7] for p > 2 and we can use this information to show 
that £{D) is bounded by powers of the volume \D\. E.g., the L°° norm is at most vr~", so 
£{D) < 7r-"|D|. 

For certain D, however, £{D) is not only finite, it is even less than 1. In [2], Flan- 
drin conjectured this to be true for all convex domains, and he showed that for all cen- 
tered two-dimensional discs B^{r) of radius r, the standard normalized Gaussian function 
Tr""*^' ■* exp(— a;^/2) is the unique maximizer of ()1.2 h . In particular £{B'^{r)) = 1 — e~'^ 
(see [2] , cf. [3] ) . It follows immediately from the definition of the Wigner distribution that 
Flandrin's proof can be easily generalized to higher dimensional poly-discs because the 
maximization problem then has a simple product structure. A less obvious case is the 2n- 
dimensional Euclidean ball B'^''^{r). The following is the generalization of Flandrin's result, 
and our main result: 

Theorem 1. The standard normalized Gaussian vr""''^ exp(— x^/2) in L2(M") is the unique 
maximizer of the Wigner distribution localization problem for any 2n- dimensional Euclidean 
ball centered at the origin. In particular, 

£ [B'-ir)) = \j e-^^M^vi) dxdy = 1 - f^, (1.3) 

where r(s,x) = J°° t^~^ e~^dt is the upper incomplete gamma function. 

Remarks: (1.) Owing to the translation covariance of the Wigner distribution, equa- 
tion (|1.3 p also applies to a ball of radius r centered anywhere in R^". It is only necessary 
to multiply the Gaussian by an appropriate linear form exp(o • x). Moreover, since the 
localization function (|1.2 |) is invariant under the action of the group of linear symplectic 
transformations, Theorem 1 can also be adapted to any image of the Euclidean ball under 
linear symplectic maps. 

(2.) Another generalization is to replace the integral over the ball with the integral 
over M^", but with a weight that is a symmetric decreasing function (i.e., a radial and 
non- increasing function of the radius ^/x^ + y^). By the "layer cake representation" [51 
sect. 1.13] the standard Gaussian again maximizes uniquely. 



2 Proof of Theorem [T] 

We start with the following preliminaries. Recall that the mixed Wigner distribution of two 
states il^i,ip2 S L^(M") is defined by 

W^i,^,(x,y) = (27r)-" / ^P^{x + T /2)r2ix - T/2)e-'^y dT . (2.1) 

Note that in contrast to (jl.l p . W^^ ^^ i^ ^^^ generally real, but, nevertheless, Hermitian 
i.e., W^^^^2 = ^^2 tt'i' Moreover, it is not hard to check that the mixed Wigner distribution 
is sesquilinear. 

Next, let /i = (^1, . . . , Hn) be a multiindex of non-negative integers, and let x G M". The 
Hermite functions Hfj_{x) on R" are defined [141 [T5] to be the product of the normalized 
one-dimensional Hermite functions, i.e., H^{x) = \Yj=i ^Mj(^j)' where 

hk{x) = n-i (kin 2-1 (-1)^ e-V2 i^e-x^ . (2.2) 

It is well known that the {i^^t} form a complete orthonormal system for L^(]R"), and that 

MH^ = \fx\H^, (2.3) 

where |/i| = X^?=iA'j, and IHI is the Schrodinger operator IHI = — ^A + ||a;p — ^. Here 
A denotes the standard n-dimensional Laplacian. In particular, the sesquilinearity of the 
Wigner distribution implies that for any ip G L2{W^), one has 

^^ = E E<^' ^m) (^. ^-)* WH„/f. . (2.4) 

fj, V 

The following lemma shows that the integral of the off-diagonal elements of (|2.4 p over any 
centered ball B'^'^ir) vanishes (cf. [5] Section 2.3). 

Lemma 2.1. Let fi, v he two multi-indices with ^i ^ v. Then, for every r ^ 0, one has 

I WH„Hjxdy = . (2.5) 

_B2"{r) 

Proof of Lemma 12. IL It is well known (see e.g. [6]) that for the one-dimensional Hermite 
functions {h^}, one has: 

fvr-i(A;!/i!)V2(_i)^^(V2zi)i-'=e-(l-l^)Li-^2|zip) ifi>fc. 

Here ^i = j;i + iyi, and L" are the Laguerre polynomials defined by 

r~"'p^ d^ 
L^ix) = — r ^(e-^-^^^), (2.7) 

for J > and a > — 1. Hence the lemma holds in the 2-dimensional case, i.e., when n = 1, 
because the integral of z^ or 1^ over any circle centered at the origin equals zero when 



j 7^ 0. The higher-dimensional case follows for the same reason from (j2.6 p . the fact that 
the Wigner distribution function yVHf^,H^{x,y) is the product of yVhm-,hni^j^yj)^ ^^^ the 
rotation invariance of the ball i?^"'(r). D 

An immediate corollary of Lemma I2.H definition ()1.2 p , and equality (|2.4 p is 
Corollary 2.2. In the notation above, 

£ (i?2"(r)) = sup / Wu.dxdy, (2.8) 

where the supremuni is taken over all multi-indices ^ = (^i, . . . ,^n) of non-negative inte- 
gers. 

The following lemma is the main ingredient in the proof of Theorem [Ij 
Lemma 2.3. For any integer A > and multi-indices fii, ^2 with A = |^i| = |/i2|, one has 

/ W// dxdy = / Wudxdy, for every r ^ . (2.9) 

Postponing the proof of Lemma 12.31 we first conclude the proof of Theorem [TJ 

Proof of Theorem [H It follows from Corollary 12.21 and Lemma 12.31 above that 

£{B^^{r))=snp f Wn.dxdy, (2.10) 

A JB^'^ir) 

where fi\ = (A,0, ... ,0), and A is a non-negative integer. Moreover, from (12.6 |) and the 
definition of the Wigner distribution it follow that: 



A 



^^H,,{x,y) = ^-^e-^'^^^^'^+y^^ L,{2ixl + yl)), (2.11) 

where L\{z) are the a = Laguerre polynomials ()2.7 p . Setting Zj = Xj -\- iyj, we conclude 
that 

WH,^dxdy= / e ^=' ' ( / ^^ e'l^^l' La(2|zi|2) dzi)dz2 • • • c?z„ 

^' E l^iP<'-2 kiP<^2-EkjP 

(2.12) 
On the other hand, from Flandrin's result in the 1-dimensional case [2], it follows that 

J Wh,dxidyi= j (-l)^e-l^il'LA(2|zi|2)dzi5; / e-l^^l' dzi, (2.13) 

S2(q,) |2i|2<a2 |2l|2<a2 

for every radius a ^ 0. An examination of Flandrin's proof reveals that the inequality is 
strict for A > 0. Hence, for every non-negative integer A one has 

WH.Jxdy^TT-'' / e ^=1 dzi---dzn = l- r^ ^ •; (2.14) 

D2n(„) n 



with equality only for A = 0. The proof of Theorem [T] now follows from (j2.ll p and (J2.14 p . 

D 

Remark: The integral in (j2.10)l is not monotone in A or in r (except for A = 0), as 
might have been thought. See [U Fig. 2] and [2] for interesting graphs of these integrals as 
a function of r. 

For the proof of Lemma [2. 31 we shall need the following preliminaries. For a non-negative 
integer A denote 

-Hx = span{Hf, ; |/i| = A} C L^{W) . (2.15) 

It follows from (12.3 p above that the space Hx consists of the eigenfunctions of the rotation 
invariant Schrodinger operator H = —2^ + fl-^l^ ~ § with eigenvalue A. In particular, it 
is a finite-dimensional, 0(n)-invariant subspace of L^(R"') with orthonormal basis {H^ : 
\fi\ = A}. It follows that for every TZ G 0{n), and every Jl with \J1\ = A, one has: 

Hji{nx)= Y. c,iil,n)H,{x), (2.16) 

u : \i'\=X 

where the coefficients Cyijl^TZ) satisfy ^ |ci,(pt,7?.)p = 1. 

We note the following useful fact: In order to identify which coefficients Cy{Ji,TV) are 
non-zero, it is only necessary to check the leading powers on the two sides of (|2.16 p . That is, 
the left side of ()2.16)l defines a polynomial of degree A in the indeterminates xi , . . . , x^. The 
highest degree terms are the monomials x'^^ ■ ■ ■ Xn" with Yl^=i f^j ~ ^^ but there are also 
monomials of degree lower than A. In order to show that a given H^, appears with a non-zero 
coefficient in the decomposition (J2.16)L it is only necessary to show that there is a highest 
degree monomial x^^ • • • x^" in the decomposition. It is not necessary to check the lower 
degree monomials; they will appear automatically because we know that the decomposition 
contains only Hermite functions of degree A and no others. 

Proof of Lemma 12. 31 Fix a non-negative integer A, and r ^ 0. We consider the maximum 
problem 



max / Wn^dxdy, (2.17) 

fJ.:\fJ.\=X J 

and denote by Jl one of its maximizers. 

From the sesquilinearity property of the Wigner distribution and Lemma 12. H we con- 
clude that for every IZ G 0{n) one has: 

j WH~{nx)dxdy = Y,\c,{jL,n)\'' j Wn.dxdy . (2.18) 

Since Hji is a maximizer, this implies that for any vq with Cin^[Ji,7l) 7^ one has 

j WH~,dxdy= j WH~(^nx)dxdy= j Wu.Jxdy , (2.19) 

_B2"{r) B2n(r) B2n(^) 



i.e., -ff^Q is also a maximizer. The lemma will be proved if we can show that, starting from 
any /x, we can, by a succession of rotations and intermediate indices, finally reach any given 
u. 

The proof will proceed in two steps. The first is to go from /I, by a succession of 
two-dimensional rotations, to (A, 0, 0, . . . , 0) with A = Yl^=i J^j- 

First, we show that there is a rotation TZ' G 0{n) with 

c^'(/i,7^') / 0) where Jl' := {{Jli +'jl2),0,Jj-3, ■ ■ ■ ijin) ■ (2.20) 

Thus, Jl' is also a maximizer. In a similar fashion, we can go from Jl' to Ji" , where Jl" := 
{{pi + /i2 + V-s), 0, 0, Ji4, . . . , Jin)- Proceeding inductively, we finally arrive at the conclusion 
that (A, 0, ... , 0) is a maximizer. 

A rotation TZ' that accomplishes the first step to Jl' is simply 7?.' : xi — )• {xi+X2)/V^, X2 — )■ 
{xi — 3;2)/\/2, Xj — )• Xj for j > 2. The monomial x'^^xi^'^ becomes ^{xi + X2)'^^(xi — ^2)^^ 
and this obviously contains the monomial xf"^ ^^ with a non-zero coefficient. 

The second step is to go in the other direction, from (A, 0, ... , 0) to (i^i, 1^2, ■ ■ ■ , '^n) when 
^^^i yj = A. As before, we do this with a sequence of two-dimensional rotations, the 
first of which takes us from (A, 0, ... , 0) to (A — 1^2, 1^2, 0, ... , 0). From thence we go to 
(A — j^2 ~ ^3) ^2i 1^3) 0, . . . , 0), and so forth. This can be accomplished with the same rotation 
as before, namely TZ' : xi ^ (xi -|- X2)/\/2, X2 -^ (xi — X2)/\/2, Xj — ;• Xj for j > 2. 

D 
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